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Abstract. We examine theoretically the behaviour of the inner accretion disk in GRS 1915+105 when soft X-ray dips 
■ are present in the X-ray light curve. We assume the presence of a radial shock in the accretion disk, as in some of the 
Two Component Advective Flow (TCAF) solutions. We discuss the behaviour of the flux tubes inside a TCAF (which 
we name Magnetized TCAF or MTCAF model for brevity) and compare various competing forces on the flux tubes. In 
^ this MTCAF model, we find that the magnetic tension is the strongest force in a hot plasma of temperature ^ 10^°K 
and as a result, magnetic flux tubes entering in this region collapse catastrophically, thereby occasionally evacuating 
\^ , the inner disk. We postulate that this magnetic 'rubber-band' effect induced evacuated disk matter produces the 
P\| blobby components of outflows and IR/radio jets. We derive the size of the post-shock region by equating the time 
scale of the Quasi-Periodic Oscillations to the infall time of accreting matter in the post-shock region and found the 
shock location to be ~ 45 — 66rg. We calculate the transition radius rtr, where the Keplerian disk deviates into a 
sub-Keplerian flow, to be ~ 320rg. Based on the derived X-ray spectral parameters, we calculate the mass of this 
' region to be ~IO^*g. We conclude that during the X-ray dips the matter in the post-shock region, which manifests 
itself as the thermal-Compton component in the X-ray spectrum, is ejected, along with some sub-Keplerian matter in 
Q ■ the pre-shock region. 

^ Key words: Accretion, accretion disks - Black hole physics - Stars: winds, outflows - Stars: individual: GRS1915+105 
. . ', - X-rays: stars 
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1. Introduction 

GRS 1915+105 has proven to be an ideal source to study in detail many of the physical concepts regarding accretion 
onto black holes. Ever since its discovery (Castro-Tirado et al. 1992), it has been continuously bright in X-rays and 
it shows a variety of X-ray variability characteristics (Morgan, Remillard and Greiner 1997; Muno et al. 1999; Yadav 
et al. 1999; Belloni et al. 2000). It has been monitored extensively in the radio band (Mirabel and Rodriguez 1994; 
Pooley and Fender 1997; Fender et al. 1999) and several episodes of high radio emissions, huge flares associated with 
superluminal motions, radio oscillations etc. are observed in this source. Several attempts have been made to associate 
the radio emission, presumably coming from jets, to the X-ray emission from the accretion disks (Fender et al. 1999; 
Naik et al. 2000; Naik and Rao 2000). The Spectral signature of winds from the Comptonising region has also been 
identified (Chakrabarti et al. 2001). 

Chakrabarti and Manickam (2000, hereafter CMOO) have applied the Two Component Advective Flow (TCAF) 
model of Chakrabarti and Titarchuk (1995) to explain various X-ray variability characteristics of GRS 1915+105. 
Recently there has been overwhelming evidence that the TCAF model is valid for many black hole candidates (Smith, 
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Heindl and Swank, 2001; Smith et al. 2001). CMOO invoked outflows from the inner accretion disk to explain a 

correlation between the QPO frequency and the duration of the burst-off states during a regular oscillations seen in 
the source. These outflows, however, are confined to the sonic sphere and fall back on the accretion disk after being 
cooled down by an inverse Compton effect. It was pointed out by Naik and Rao (2000) that the source does not show 
appreciable radio emission during such oscillations. A detailed inflow/outflow model has not been presented for this 
source to explain the radio emission, particularly to explain the superluminally moving radio blobs. 

Recently Naik et al. (2000) have detected a series of soft X-ray dips during the declining phase of a huge radio flare 
and have postulated that such soft dips are responsible for the jet emission. Vadawale et al. (2001) made a detailed 
study of X-ray dips observed during the radio flare using the Rossi X-ray Timing Experiment (RXTE) data and have 
presented evidence for the disappearance of the inner accretion disk during the dips. Since the disappearance of the 
inner disk is seen to be correlated with intense radio activity, the role of the magnetic field must be studied in order 
to understand the system completely. Rodriguez and Mirabel (1999) estimated the field in radio blobs to be around 
tens of mG at 500-lOOOAU (in 1994 observation). Fender et al. (1997) requires the field to be around 8G at around 
lAU (in their 1996 observations). From the similarity of ^ 30min oscillations in IR and Radio, they concluded that 
the radio blobs are adiabatically expanding and are independently ejected from the disk every 30 — 40 minutes. If the 
trapped fleld inside a radio blob is of roughly 1 jr (for toroidal field) then its interpolated value close to a black hole is 
around lO^G at around lOr^ which is comparable to an equipartition value. Thus, one needs to correlate fields ejected 
from the disk with those observed inside the radio blobs. In this Paper, we examine the mass ejection based on the 
TCAF model in presence of a magnetic fleld (we call this as Magnetized TCAF or MTCAF model) amplified due to 
strong shear at the transition radius of the Keplerian and a sub-Keplerian flow. In the next Section, we discuss forces 
which govern the motions of the flux tubes and write equations of toroidal flux tubes inside an accretion disk with a 
constant angular momentum. We show that close to the black hole, where the flow could be very hot ( ^ 10^°K)the 
flux tubes move at least with the Alfven speed and may catastrophically collapse like a stretched rubber band. We 
conjecture that such rapid collapse would assist evacuation of matter from the disk and cause X-ray 'dips' seen in 
the light curves. In §3 we estimate the mass of the ejecta which agrees with observations. Finally, in §4 we draw our 
conclusions. 

2. The Magnetized TCAF model for GRS 1915+105 

Based on the global solutions of the most general advective accretion disk solution (Chakrabarti, 1990, 1996a), 
Chakrabarti and Titarchuk (1995) presented a TCAF model of accretion onto black holes. According to this model, 
matter with high viscosity and angular momentum settles into the equatorial plane, while matter with lower viscosity 
and angular momentum surrounds the Keplerian disk. This sub-Keplerian component is likely to form a standing or 
an oscillating shock (Molteni, Sponholz and Chakrabarti, 1996; Ryu, Chakrabarti and Molteni, 1997) front near the 
centrifugal barrier (few tens of Schwarzschild radii) depending on whether the Rankine-Hugoniot condition is satisfied 
or the cooling time in the post-shock region is comparable with the infall time. Soft photons from the Keplerian 
disk in the pre-shock flow are intercepted by the puffed up sub-Keplerian post-shock flow and are reprocessed due 
to inverse Comptonization. If the post-shock matter remains hot, the black hole exhibits a harder spectrum, while if 
the post-shock region is cooled down by the soft photons, the black hole exhibits a soft spectrum. There could be a 
sub-Keplerian region just before the shock as well since with low viscosity and accretion rates a Keplerian disk recedes 
from a black hole and it is not necessary that the shocks form right at the region where a transition from Keplerian 
to sub-Keplerian flow takes place. Fig. 1 schematically shows this behaviour of the two components. 

The centrifugal pressure supported boundary layer (or, CENBOL for short) formed in a transonic, advective flow 
(Chakrabarti et al. 1996), has most of the features of a thick accretion disk although in advective disks, advection 
is included self-consistently and in thick disks advection is totally ignored. This is because at the CENBOL surface, 
matter undergoes a supersonic to sub-sonic transition and it moves very slowly in the radial direction. In some phases 
of accretion, matter can bring in a large stochastic magnetic field. The field is sheared due to strong azimuthal velocity 
and the toroidal field becomes very strong. These field lines will have very little matter within it and would likely 
be buoyant and emerge from various parts of the accretion disk. Fig. 1 schematically shows this behaviour. Toroidal 
magnetic flux tubes released from the Keplerian disk are sheared, amplified and are advected in the sub-Keplerian 
fiow. Due to the centrifugal barrier matter stays away from the axis. Thus a so-called funnel wall is created exactly 
as in a thick accretion disk (Paczyhski & Wiita, 1980). Chakrabarti and D'Silva (1994, hereafter CD94) computed 
the nature of their trajectories inside a thick accretion disk. They showed that in the event a strong flux tube enters 
a hot region with ion temperature Tj ^ lO^^K, the magnetic tension becomes the strongest force and the flux tube 
catastrophically collapses. Because of strong similarity of the thick accretion disk and an advective flow, especially 
inside the CENBOL, we believe that a similar mechanism could be working and flux tube collapse would take place. 
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We shall consider the motion of the flux tubes on the equatorial plane of an accretion flow around a Schwarzschild 

black hole described by Paczynski-Wiita (1980) pscudo- Newtonian acceleration g = — 1/2(.t — 1)"^. Wc use the 
geometric units. Masses are measured in units of the mass of the central black hole, Mbh] distances from the axis 
{x) will be measured in units of the Schwarzschild radius rg = 2GMbh/(?', and the time scales are measured in 
units of Tg/c. Inside the disk, we choose the polytropic equation of state. P = Kp^ , where K and 7 are constants. 
The magnetic flux tubes brought in by advection are assumed to be sheared and axisymmetric toroidal flux tubes of 
random shape and size could be produced inside the flow. We however assume that the flux tubes are thin, i.e., the 
flux tube cross-sectional radius cr is smaller than the local pressure scale height of the disk. Close to a black hole, 
angular momentum of the flow remains constant (Chakrabarti, 1996a) even in presence of moderate viscosity. Thus, 
we choose speciflc angular momentum A in the sub-Keplerian region to be constant. The equations of motion for thin 
flux tubes have been given in CD94 and we do not repeat here. For the sake of completeness, we write down the radial 
equation only valid for the equatorial plane {9 = 7r/2), 

where X — rm/me, rm — 27r^a'^xpi and m,. = 2'K^a'^xpe are the masses of the fluid inside and the fluid displaced 
by the flux tube respectively, pi and pe being the corresponding densities. Subscripts e and i indicate whether the 
relevant quantity is of the ambient (external) medium or within the flux tube (internal medium), tp = wa^B, B being 
magnetic fleld of the tube and 4> is the intrinsic angular velocity of the flux tube inside the disk. The drag term is 
assumed to be similar to the drag experienced by a cylinder moving perpendicular to its axis in a fluid, and the drag 
per unit length in radial direction is given by, 

D^ = ~CDpe(y{x-uf (2) 
where Cd ^ 0.4 (Goldstein 1938). The pre-shock, sub-Keplerian flow is assumed to have a radial velocity. 



M~/3/V(x-l) (3a) 

and the post-shock sub-Keplerian flow, 

u~ l/i?v'(a;- 1). (36) 

Here, /? is a factor by which sub-Keplerian matter slows down compared to a freely falling flow. /3 ^ 1 for very low 
angular momentum cool flow. R is the compression ratio of the shock by which matter is assumed to be slowed down 
inside the CENBOL. 

Prom Eq. (1), we note that there are four forces in operation: (i) The second term inside the bracket of the left 
hand side is the Coriolis force Fc = 2v^ui, where u = X/x^ is the angular velocity of the flow. The Coriolis acceleration 
is given by (CD94), 

Pi 2A2 , 1 1 

ac = — [— 2 J- (4) 

Pi \ Pe X Xq X 

A flux tube brought from xq to x inside a flow of constant angular momentum would feel no Coriolis force if the 
motion is along the direction of constant angular velocity since ~ a; for rotating bodies. A magnetic flux tube with 
buoyancy factor Mi,{= P''~P' ) = 1 — X, will feel the (ii) Magnetic Buoyancy force (flrst term on the right hand side 
inside the curly bracket). The corresponding acceleration is (CD94): 

aMB-Y^[2(x-l)2-^]- 

The tube will emerge out of the disk if cimb dominates over ac- The middle term on the right hand side of Eq. (1) is 
the (iii) force due to Magnetic Tension. The corresponding acceleration is (CD94): 

1 _ 
rui + rUe 2iTa'^ A'Kx{pi + pe) 

The flnal and very important force is represented by the final term of Eq. (1). It is called the (iv) Drag Force exerted 
due to the motion of a rigid tube inside a fiow. The corresponding acceleration is (CD94): 

1 CDPia-{x - uf CnPiV^ 

Oj d . ( I ) 

2 27ro-2(pi + pe) 2na{pi+pe) ^' 



4 



Nandi et al: Ejection of inner accretion disk in GRS 1915+105 



Here we used i for computing the maximum value of the drag force. 

Equating the acceleration due to buoyancy with that due to magnetic tension (Eq. 6) on the equatorial plane, we 
get the critical temperature of the external gas as (CD94) , 

TOpC^[A|^(xo)-A^] 

above which the tension dominates over buoyancy and the flux tubes rush towards the funnel wall catastrophically. 

Here we used the mean electron number per ion to be /i = 1/2 and k is the Boltzmann constant. Subscript specifically 
indicates that the flux tube will behave like an over stretched rubber-band only after it crosses x = xq, where T > Tp o- 
Note that since we are dealing with a sub-Keplerian disk, Tp^ > in the entire region of interest. For typical values 
xq = 100, Ao = 1.8, one obtains Tp^o ~ 6 x lO^K. For an adiabatic disk, this assumption remains valid even when the 
flux tube collapses very rapidly since its internal temperature will increase adiabatically in the same way as in the 
external disk. For an efficiently cooled two-temperature flow, the above proton temperature would correspond to an 
electron temperature of Tg ^ y/ nie/mpTp. This is around 14keV which is very reasonable for the temperature of the 
sub-Keplerian region. When the flux tubes fall radially, neither Coriolis force nor the drag term could be neglected. 
In fact, accelerated flux tubes would have high x in the drag term as they move faster than the bulk radial motion. 
Meanwhile, assuming that internal mass of a flux tube is roughly constant, the cross-section a ^ ^l\fx increases as 
the flux tube approaches the black hole. So, this will also increase the drag term. As a result, we expect that the 
flux tube would slow down somewhere close to the black hole and buoyancy would eject the flux tube out of the disk 
perpendicularly as shown in Fig. 1. Typical trajectories of the flux tubes, based on numerical integrations in CD94, 
are shown in Fig. 1. 




Fig. 1. A cartoon diagram of the accretion disk near a black hole which includes a shock (rs), a sub-Keplerian and a Keplerian 
disk with a boundary at rtr. Stochastic magnetic fields are sheared and amplified as they leave a Keplerian disk. In a hot, 
sub-Keplerian flow, these toroidal flux tubes catastrophically collapse squirting matter along the axis, and thereby evacuating 
the disk and producing outflows. Dashed curves show typical trajectories (CD94). 

The buoyancy timescale 4 is estimated from Eq. (5): auB ~ ^/^^ = {\^x)x'^ ~ where, h ~ agX'^/^ is the 
half thickness of the disk in vertical equilibrium at x and ag ^ x~^l'^ is the speed of sound. Note that t\, is very large 
in a Keplerian disk {\k ^ A). For a sub-Keplerian flow, A/Ak 0.4 with A = 1.8 at a; = 30, X 0.1, 0.9, 
U ~ 240. 

Collapse time of a flux tube is estimated in the following way: When tension is the most dominant force, the 
radial equation of the flux tube (Eq. 1) is simplified to: 

^ -I- At = 0, (9) 
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Putting X x/t'j, we see that the velocity of collapse of the flux tube Vf is, 

where Va is the Alfven velocity. Thus, flux tube collapses in Alfven speed. Since B oc x~^, magnetic pressure is oc 
while the gas pressure oc x^^^^ , the ratio of magnetic to gas pressure / cx x^^^^ , increasing with the decrease of x. 
As the flux tube leaves the Keplerian disk, large shear at the transition region and in the sub-Keplerian flow quickly 
intensifies the flux tube to M(, ~ 1 during infall. On the other hand, if / ~ 1 at x ~ 300 where the flow deviates 
from a Keplerian disk, / ~ 30 at a; = 30, inside the CENBOL. The Alfven speed Va ~ (jf^)^/^ ~ if/lV^^as ~ Sa^ 
where we used 7 = 4/3. Since, inside a CENBOL, the velocity of matter is highly sub-sonic, it is therefore high 
sub-Alfvenic: Vm « Va- Thus, Vf/vm » 1 and the flux tube collapses catastrophically. This justifies the conjecture 
made earlier that the rubber-band effect could evacuate the disk (Chakrabarti, 1994; 1996b; 2000). For ~ x~^^'^, 
tf ~ x3/2/5 ~ 30 at x = 30. Thus inside the CENBOL tb » tf. 

So far, we have ignored many non-axisymmetric effects such as Parker instability (1979 and references therein) and 
shearing instability (e.g., Balbus & Hawley, 1991). Foglizzo & Tagger (1995) treated this problem comprehensively in 
the context of a standard disk embedded in a large scale field. They found that (a) if the wave-length is larger than the 
disk thickness then the flux tube is very unstable and buoyantly comes out of the disk and (b) instability is strongest 
if the fleld is weaker. According to Parker (1979), submerged field tubes may break up into filaments in timcscalcs 
of around 3A/va, where A is the scale height and the field is able to escape from the gas in timescales tp A/va- If 
A ^ /i ~ a;, which is especially true in CENBOL region, the time scale of the escape of the field may be comparable to 
the tf as obtained above. As a result, the flux not only collapses catastrophically, but also escapes upwards following 
a curved trajectory as depicted in Fig. 1. In presence of a differential rotation, Balbus and Hawley (1992) suggested 
that even a small initially vertical fleld would be amplified to create all the components in the dynamical timescale 
^ l/w = x^/X ~ 50 (at X = 30, A = 1.8). Numerical simulation (Hawley, Gammie & Balbus, 1995) has verified 
this instability neglecting the tension effects. Even though the time scales of the collapse, the Parker instability and 
the shear instability are of the same order, we believe that the rubber-band effect would still be important for the 
destruction of the inner disk. 



3. Estimation of mass of the post-shock region and the sub-Keplerian region 

Vadawale et al. (2001) have shown that during the soft X-ray dips a thermal-Compton component in the X-ray 
spectrum gets suppressed. Several works in the literature talk about the disk-evacuation (e.g., Belloni et al. 1997; 
Feroci et al. 1999) in this context. We like to imdcrstand this using Two Component Advcctivc Flow (TCAF) models 
of Chakrabarti & Titarchuk (1995) in presence of a Magnetic field (i.e., MTCAF model) and its time variability 
properties described in Molteni, Sponholz & Chakrabarti (1996); Ryu, Chakrabarti & Molteni (1997) and Chakrabarti 
& Manickam (2000). The observation of possible disk-evacuation is clearly in line with the TCAF model and the 
shock oscillation model of the quasi-periodic oscillations (CMOO; see also, Rao et al. 2000) which showed that the 
Comptonising post-shock region participates in oscillation. The sub-Keplerian region in the pre-shock flow does not 
emit miich radiation and it is possible that some of this region may also be disrupted during the rapid collapse of the 
flux tube. Once we accept the destruction of the sub-Keplerian region by the magnetic rubber-band effect, we can 
compute the mass of this region in the following way: 

The shock location is computed by equating the infall time from the shock with the time scale of QPO. This time 
scale tff can be written in the form (CMOO): 

tn = i-—^ (11) 
Rr^ rg' ^ ' 

where, R is the compression ratio (see also, Eq. 3b), is the shock location, vq is a dimensionless constant. Here, 
a = 3/2 for free-fall motion and a = 1 for a flow of constant velocity cvq/R in the post-shock region. Using this 
assumption, the shock location in the pre-dip and the post-dip flow (for parameters in these flow see, Vadawale et al. 
2001) are ^ 45rg and ^ 66rg respectively. CMOO proposes that a better fit in the correlation between the diuation 
of the QPO and the frequency of QPO requires more or less constant velocity in the post-shock region with a rough 
velocity of 0.066c/i?. This produces the shock in the pre-dip and the post-dip flow at 79rg and 140rg respectively. 
Using the parameters of Vadawale et al. (2001), and the location of the shock as given above, the electron number 
density in the Comptonising region is found to be Ug ~ 10^^ cm~^. The corresponding mass {^nr^nipnerg) of the 
region is Mcenbol ^ 2 — 3 x lO^^g depending on the model of the inflow. How much of this matter is squirted out of 
the disk along the axis? Again, in the absence of the size distribution of the flux tubes, the answer is difficult. However, 
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because of centrifugal force, matter is unlikely to enter within the 'funnel wall' (CD94) even after compression, sudden 
collapse of large flux tubes [a ~ h{rs)] would be expected to displace the whole CENBOL region parallel to the funnel 
wall. The outflow rate would be Mout ~ McENBOh/tf which may be very large compared to the inflow rate for a 
short duration of t/. If flux tubes are smaller in size, since the collapse velocity is much larger than the speed of 
sound, matter will still be displaced but would be refilled in free-fall time, unless there are many near-simultaneous 
flux collapse events. 

If the magnetized sub-Keplerian disk is removed by imploding flux tubes as described in the earlier Section, one 

requires to know the location of the inner edge of the Kcplcrian disk to estimate the complete mass involved. From 
the model fit (Vadawale et al. 2001), the Keplerian disk temperature turns out to be Tk ~ 1.5keV. With a hardening 
factor of around 1.7 (Shimura and Takahara, 1995), the mass of the black hole as lOM© and Shakura-Sunyaev 
viscosity parameter 0:55 = 0.01, the above temperature corresponds to a transition radius at around rir ^ 320rg 
(Shakura and Sunyaev, 1973). Assuming density falling off as p ~ po(^)^^^) the mass of the sub-Keplerian flow of 
size xtr (/ 4:TTxh{x)pdx with h{x) x) is around 10^'^g. These computations assume no pair production, i.e., there is 
exactly one electron for each proton in the Comptonising region. 

Once the evacuation is complete, the disk is filled in quickly by sub-Keplerian matter in timescale of: 

t.s. - ^i^y^^ = 192(^)(^)^(||-)V^-|^. 
ass X VKep oiss a-s 300 lOM© 

Here, h{x) ^ agX^/"^ {as is the speed of sound) is the local vertical height of the disk. The time scales seem to be 
reasonable, since the dips have been seen to be filled up in a matter of 150 to 200 seconds. It is to be noted that 
Dhawan et al. (2000) using disk-instability model of Belloni et al. (1997) obtained the missing inner disk region to be 
only 180km. This distance is about 6rg and with lO^^gm s"-'^ accretion rate which they employ, the mass of this region 
cannot be enough to create the lO'^'^gm blob (or, even lO^^g mini-blobs) observed by Mirabel & Rodriguez (1999). 
This is particularly because the region 1 — 3rg is definitely supersonic and sub-Keplerian and therefore is of very little 
mass. We believe that the missing region should be much larger, possibly order of a 100 rg or so. 

4. Conclusion 

In this Paper, we have given a physical basis for a sudden mass ejection in GRS 1915-1-105. We showed that if 
matter brings in a particularly strong magnetic field, this would be sheared and amplified to a value much above the 
equipartition value before it can be expelled by buoyancy. Magnetic tension collapses these toroidal flux tubes at a 
highly supersonic speed, much faster than the flow velocity. This has the effect of displacing matter from the disk in 
transverse direction (much like a fast boat causing spillage on a shore) and depositing it to outflowing winds. From 
the observed fits of Vadawale et al. (2001) we estimated the electron number density and the mass of the post-shock 
region and the sub-Keplerian flow to be around lO^^g and 10^°g respectively. Our estimate of the post-shock mass is 
a factor of ten less than the mass estimate (lO^^g) of 'baby-jets' (Mirabel et al 1998) associated with IR and radio 
flares and could therefore be due to ejection of some sub-Keplerian matter as well. These 'baby-jets' are found to be 
associated with class /3 light curves which have soft X-ray dips. These dips are also seen in class 9 light curves. During 
a major portion of the huge radio flares associated with superluminal blob emission a series of soft dips are present 
(Naik et al. 2000). Mirabel and Rodriguez (1999) have pointed out that in each epoch of this type of outflow, the mass 
condensation is around lO^^g. In order to achieve this, we require that matter is accumulated from disk evacuation at 
each 'dip' and within each epoch, successive mini-blobs move faster than the earlier blob in order to 'catch up'. This 
may indicate some other runaway process with a longer time scale of tens of days. Naik et al (2000) have observed 
such X-ray dips at a rate of once in a few hundred seconds during the peak or the radio flare. If there are ~ 1000 
evacuation events during an episode of superluminal blob ejection (in a few days), then total mass condensation would 
be lO^^g. Hence in order to explain the observation of Mirabel and Rodriguez (1999), one must reqiure that in each 
epoch, matter is accumulated from at least a thousand evacuation events. Future observation would tell if such is the 
case. 

It is to be noted that the mass of the condensation as estimated by Mirabel and Rodriguez (1999) is based on the 
presence of one electron per proton, i.e., no pair production is assumed. With a pair density, say, ten times larger, the 
estimated mass would be ten times less. However, at the same time, estimated mass of the disk would also be reduced 
by the same factor. Hence, the number of ejection events is not affected. 
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Astrophysical Flows Around Compact Objects. The authors thank the referee for helpful comments. 
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